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Conclusions: On the basis of our initial analysis, we conclude 

that the observed pattern of volcanic features may be correlated with 
the distribution pattern of global physiographic and geologic char- 
acteristics. The distribution of volcanic centers and regional tec- 
tonic patterns suggests that volcanic features are generally excluded 
from lowlands and regions of tectonic shortening, and occur pre- 
dominantly in upland regions characterized by geologic evidence 
for extension. Three hypotheses that may account for the observed 
distribution and geologic association may be categorized as 

(1) environment/elevation-related, (2) mantle dynamics -related, 
and (3) age-related. It is likely that all three influences occur, but on 
the basis of the global association of areas of high volcanic center 
abundance with tectonic characteristics of extension and the prob- 
able association of many individual volcanic centers with local 
mantle upwelling and plumes, we believe that the regional concen- 
trations of volcanic centers may be primarily associated with 
regions of broad mantle upwelling phenomena. Although the broad- 
scale characteristics and association of the distribution of volcanic 
centers may be accounted for by the first hypothesis, details of the 
distribution and local associations may be strongly influenced by 
altitude and age-dependent effects. 
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THE SPIN VECTOR OF VENUS DETERMINED FROM 
MAGELLANDATA. M. E.Davies 5 ,T. R. Colvin 1 , P.G. Rogers 1 , 
P. W. Chodas 2 , and W. L. Sjogren 2 , 3 RAND, USA, 2 Jet Propulsion 
Laboratory, Califomialnstitute of Technology, PasadenaCA91 109, 
USA. 

A control network of the north polar region of Venus has been 
established by selecting and measuring control points on full- 
resolution radar strips. The measurements were incorporated into a 
least-squares adjustment program that improved initial estimates of 
the coordinates of the control points , pole direction, and rotation rate 
of Venus. The current dataset contains 4206 measurements of 606 
points on 619 radar strips. The accuracy of the determination is 
driven by spacecraft ephemeris errors. One method used to remove 
ephemeris errors is to adjust the averaged orbital inclination and 
argument of periapsis for each orbit. A more accurate method that 
has been used with selected blocks of orbits incorporates optimally 
fitting measurements of additional points at all latitudes of Lhe radar 
strips together with Earth-based spacecraft radiometric tracking 
measurements to compute new spacecraft ephemerides. The root- 


mean-space (RMS) of the point measurement residuals in these 
improved ephemeris solutions is typically about 20 m in slant range, 
and 40 m in the along-track direction. Both the control network 
computations and the improved ephemeris solutions incorporate 
radii at the measured points derived from the Magellan altimetry 
dataset [1]. The radii of points north of 85° are computed in the least- 
squares adjustments. 

An accurate estimate of the rotation period of Venus was 
obtained by applying the ephemeris improvement technique to the 
second cycle closure orbits 2166—217 1 that overlaid the first cycle 
initial orbits 376-384. Sixty-four common points were measured on 
both orbit groups and improved ephemeris solutions computed over 
both blocks simultaneously, along with the rotation Tate. A similar 
analysis was made using orbits 874—878 from cycle 1 and 4456-4458 
from cycle 3 . Fifty-two common points were measured on both orbit 
groups and the rotation period of 243.0185 ± 0.0001 was computed. 
This latter solution confirmed the initial solution, and was an 
improvement over the first closure solution because of the longer 
period between overlapping orbits. 

The geodetic control network uses measurements of points on 
overlapping radar strips that cover the north polar region; these are 
only the even-numbered orbits. These strips were taken in the first 
cycle and encircle the pole except for three gaps due to the superior 
conjunction data loss, the reduced data due to occultation, and the 
area of ongoing work. Improved ephemeris solutions for 40 orbits 
(376-384, 520-528, 588-592, 658-668, 1002-1010, 1408-1412, 
1746-1764, and 2166-2170) are included and fixed in the geodetic 
control computations, thus tying the network to the J2000 coordi- 
nate system. The argument of periapsis and orbital inclination of all 
remaining orbits were allowed to vary as part of the least-squares 
adjustment. The RMS of the point measurements is typically on the 
order of 75 m in both along-track and cross-track. The rotation 
period was fixed at 243.0185 days. The coordinates of the 606 
measured points were determined and the solution for Lhe direction 
of the north pole was a = 272.76° ± 0.02°, 6 = 67. 1 6° ± 0.0 1° (J2000). 
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MONTE CARLO COMPUTER SIMULATIONS OF VENUS 
EQUILIBRIUM AND GLOBAL RESURFACING MODELS. 
D. D. Dawson 1 . R. G. Strom 1 , and G. G. Schaber 2 , •University of 
Arizona, Tucson AZ 85721 , USA, 2 U.S. Geological Survey, Flag- 


staff AZ 86001 , USA. 


Two models have been proposed for the resurfacing history of 
Venus: (1) equilibrium resurfacing and (2) global resurfacing. The 
equilibrium model [1] consists of two cases: In easel areas <0.03% 
of the planet are spatially randomly resurfaced at intervals of 
<150,000 yr to produce the observed spatially random distribution 
of impact craters and average surface age of about 500 m.y., and in 
case 2 areas >10% of the planet are resurfaced at intervals of >50 
m.y. The global resurfacing model [2] proposes that the entire planet 
was resurfaced about 500 m.y. ago, destroying the preexisting CTater 
population and followed by significantly reduced volcanism and 
tectonism. The present crater population has accumulated since 
then, with only 4% of the observed craters having been embayed by 
more recent lavas. 

To test the equilibrium resurfacing model we have run several 
Monte Carlo computer simulations for the two proposed cases. For 
case 1 we used a constant resurfacing area of 0.03% of the planet 
with a constant thickness and a constant 150,000-yr time interval 
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Fig. 1 . (a) Observed Venus impact crater distribution, (b) 0.03% planet equilibrium resurfacing Monte Carlo simulation, (c) Observed 
volcano size distribution equilibrium resurfacing Monte Carlo simulation, (d) 10% equilibrium resurfacing Monte Carlo simulation. The 
unfilled circles are pristine craters, the solid circles arc flooded craters, and the solid triangles are possibly flooded craters. Flooded craters 
not shown in (d). 


between resurfacing events. Another more realistic simulation used 
the observed size distribution of volcanic features determined by 
Head et al. [3] and a variable lime interval between resurfacing 
events that depended on the size of the event. Each simulation was 
run for a simulation period of 3 b.y. We took into account ejecta 
blankets by doubling the size of the craters. The center of each crater 
was marked with a dot, and the amount of the crater covered by 
resurfacing events was determined. Craters that were entirely 
covered by a single resurfacing event were considered destroyed. 
The minimum amount of flooded craters were those craters that had 
ei thernot had their central point covered or had not had an area equal 
to 100% of their area covered. The maximum amount of flooded 
craters included the minimum amount and those craters that ful- 
filled both of the above criteria but had not been entirely covered by 
a single event. Some fraction of these latter craters may have been 
destroyed, but the others were probably intact and flooded. In fact, 
this category of crater only constituted about 15% of the total 
number of craters remaining after the simulations. Furthermore, a 
constant thickness resurfacing event is unrealistic and will produce 
more destroyed craters than a more realistic tapered event that will 
produced more partially flooded craters. For each simulation the 
input crater size- frequency distribution was the observed differen- 
tial -0.2 slope for craters between 2 and 30 km diameter and the 


observed differential -3 slope at diameters greater than 30 km. After 
the simulations, however, the size distribution had changed to a -2 
slope at diameters greater than 30 km, because over the simulation 
lime of 3 b.y. the smaller craters were preferentially erased relative 
to the large craters to produce the -2 slope. Therefore, it was 
necessary to input an artificial size-frequency distribution with a —4 
slope for craters greater than 30 km diameter in order to arrive at a 
-3 slope at the end of the simulation. 

For the case 1 equilibrium resurfacing model to be valid the 
resulting crater distribution should be spatially random and there 
should be only about 4% partially flooded craters. Figures lb and lc 
show the results of the simulations compared to the observed crater 
distribution (Fig. la). Although the crater spatial distribution ap- 
pears to be random, the amount of flooded craters is a minimum of 
24% and maximum of 40% in the 0.03% resurfacing model, and a 
minimum of 30% and maximum of 46% in the more realistic 
observed volcano size distribution model. This is between 6 and 1 0 
times the amount of flooded craters observed. The computer simu- 
lations do not take into account the more than 300 splotches that are 
part of the impact record and are not embayed by lava. This would 
only exacerbate the problem by producing many more flooded 
impact events. Furthermore, the input size frequency distribution 
above 30 km diameter required to produce the observed distribution 
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is completely unlike any size distribution on the terrestrial planets. 
Therefore, the computer simulations strongly -toko* that, hecase 
1 equilibrium resurfacing model is not a valid explartnt.ori for 
resurfacing history of Venus. The observed nowand^dKOTb^on 

(210% resurfacing 

areas) simply wiU not work, except at the 100% (global) resurfacing 
level Figure Id is a Monte Carlo simulation for constant 10% 
Si! c£g areas with a constant 50-m.y. time interval between 
events. Clearly the crater distribution is completely "®» dom “ d 
totally different from the observed distnbut.on. We have done 
simulations for 25% and 50% resurfacing areas with simll = Tre ^ , 
Therefore, the equilibrium resurfacing model is not a valid model 
for an explanation of the observed crater populauon charactensti 
or Venus’ resurfacing history. 

The global resurfacing mcxlel is the most likely explanation for 
the characteristics of Venus' cratering record. The amountof 
resurfacing since thatevent, some 500 m.y. ago. can be estimated by 
a different type of Monte Carlo simulation. In this simulation die 
cratering r<2rf begins 500 m.y. ago with the observed crater size 
distribution. Our first simulation randomly selects craters from tins 
si'zedisfribulion and resurfaces areas with volcanos randomly 
selected from the observed volcano size distribution. The tun 
interval between volcanic events is chosen so that only 0 0 * 

craters are flooded at the end of 500 m.y. To date, our iniual 
simulation has only considered the easiest case to implement. In this 
case the volcanic events are randomly d.stnbuted 
planet and. therefore, contrary to observation, the flooded craters 
Ire also randomly distributed across the P 1 ™*' 
results in a maximum resurfaced area of about 10% of die planet 
since die global event, and an obliteration of about 4% of the craters^ 
Future simulations will take into account the observed nonrandom 
distribution of Hooded craters and , therefore . ^ 
bution of volcanic events. These simulations will probably result m 
a lower percentage of planet resurfacing because volcamsm will b 

“reSJsVTrSs R. J. et al. (1992) LPSC XX///, 
1065-1066. (21 Schaber G. G. et al. (1992) JGR. special Mage an 
issue, in press. [3] Head J. W. et al. (1992) JGR , special Magellan 
issue, in press. ^ r 
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MF THANE measurement by the pioneer VENUS 
large probe neutral mass spectrometer, t m 

^nahuei ^ R. R. Hodges *.». -University of Michigan. Am, 
Arbor MI 48109, USA, ^University of Texas at Dallas, Richardson 
TX, 75083, USA. 

The Pioneer Venus Large Probe Mass Spectrometer detected a 
large quantity of methane as it descended below 20 km in the 
atmosphere of Venus. Terrestrial methane and'*Xe. ' 
ing in the same container and flowing through the same plumbing, 
were deliberately released inside the mass spectrometer for instru- 
mental reasons. However, the >*Xe did not exhibit behavior sun, ar 
to methane during Venus entry, nor did CH 4 in laboratory Simula^ 
tions The CH 4 was deuterium poor compared to Venus wa er and 
hydrogen. While the inlet to the mass spectrometer was clogge 

with sulfuric acid droplets, significant deuteranon of CH.and its H 2 

progeny was observed. Since the only source of deuterium identdi- 
able was water from sulfuric acid, wehave concluded that we shoul 


correct the HDO/H 2 0 ratio in Venus water from 3.2 x 1<H to 

(5 "wten^the probe was in the lower atmosphere, transfer of 
deuterium from Venus HDO and HD to CH 4 can account quanuu- 
Cively for the deficiencies recorded in HDO and HD below . 
and consequently, the mysterious gradients in water vapor an 
hydrogen mixing ratios we have reported. The revision in the DW 
ratio reduces the mixing ratio of water vapor (and Hj) reported 

^wT^nmyeTable to say whether the methane delected was 
atmospheric or an instrumental artifact. If it was atmospheric ns 
release must have been episodic and highly localized. Otherwise^ 
die large D/H ratio in Venus water and hydrogen could not be 

maintained. a n 
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VISCOELASTIC RELAXATION OF VENUSIAN CORONAE 
S MOUNTAIN BELTS: CONSTRAINTS °N GLOBAL 
HEAT FLOW AND TECTONISM. I. Duncan and A. Leith. 
Department of Earth and Planetary Science, Washington Umver- 
sity, St. Louis MO 63130, USA. 

Venus differs from Earth in that water is essentially absent and 
its surface temperatures are about 470 K higher. Hie competing 
effects of high surface temperature and dry lithologies on die long- 
term history of surface topography have been studied using the 

finite-element method (Tccton) Ul- 

The relaxation history of surface topographic features, such as 
coronae and mountain belts, is a function of thermal gradient, 
crustal thickness and lithology . regional stresses and basal tracuons 
applied to the lithosphere. In this study we haveexammed the 
relative effects of these factors over a period of 500 M a (presume 
to be the mean age of the venusian surface) [2] 

We assume that the venusian emst is composed 
combinations of diabase, gabbro. komatiite. and refractory ,thoh>- 
gics such as anorthosite and websterite. Using appropnate thermal 

— «•» «*“ E “T 

(with and without a secular cooling contribution from * c “ re ' 

[3,4], thermal gradients ranging from about 20 K km to 

are computed. We further assume that the thickness of a diabase 
crust is limited by the dry solidus. The models are dynamically 
isostatically balanced, using an elastic foundation. 

Preliminary results of the study are shown in Fig. 1 , m which a 
2-km-high volcanic plateau has been instantaneously emplaced on 
Uiesurface. For this model the crustal thermal gradient was 
38 K km-' After the elastic response (essentially representing 
initial isostqtic balance) the topography relaxes until 
about 230 m above the surrounding region, and the slope tomn _g 
crest to moat has been reduced from an initial 6 to about 2. • 

values we obtain for our model plateau heights and slopes are in the 
observed range for venusian coronae. Thus we argue that coronae on 
Venus can be modeled as the product of elastoviscous relaxation of 
volcanic plateau. Although our starting models are oversunp ifica- 
tions they do show all the critical morphological characteristics of 
venusian coronae. Matching die observed spectrum of coral a 
morphology by vary ing the size, initial slope, and rheology of model 
plateaus enables constraints to be placed on plausible values for h 
venusian heat flux. We argue that the mean global heat flux mustbe 
significantly lower (-»« to be consistent with the observed spec- 

trum of coronae topography. , p •« j 

We are presently examining models similar to those 
above to investigate venusian mountain belts. Our models differ 



